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"^ This paper describes a general program for deriving the action of single field inflation 

\^ models with nonstandard kinetic energy terms using CMBR power spectrum data. This 

^^ method assumes that an action depends on a set of undetermined functions, each of 

which is a function of either the inflaton wave function or its time derivative. The 
scalar, tensor and non-gaussianity of the curvature perturbation spectrum are used to 
derive a set of reconstruction equations whose solution set can specify up to three of the 
undetermined functions. The method is then used to find the undetermined functions 
in various types of actions assuming power law type scalar and tensor spectra. In 
^, actions that contain only two unknown functions, the third reconstruction equation 

implies a consistency relation between the non-gaussianty, sound speed and slow roll 
^^ parameters. In particular we focus on reconstructing a generalized DBI action with 

(^ an unknown potential and warp factor. We find that for realistic scalar and tensor 

CN spectra, the reconstructed warp factor and potential are very similar to the theoretically 

derived result. Furthermore, physical consistency of the reconstructed warp factor and 
potential imposes strict constraints on the scalar and tensor spectral indices. 
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\-{ Since the landmark COBE experiment, the study of the cosmos has entered a new age of 
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precision cosmology. For the first time in the history of modern cosmology, direct quantitative 
measurements of early cosmological observables were available. The data taken by COBE 
was critical in establishing inflation as the central paradigm in our theories of the origin of 
the universe [Ij . Thanks to experiments that measured the spectrum of CMBR fluctuations, 
we have now conflrmed that the near-scale invariance of large scale fluctuations that is a 
prediction of inflation are in fact borne out in the data. Although observation supports 
the general theory of inflation, as of now the data is unable to conclusively determine the 
mechanism responsible for inflation. 
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The difficultly in discriminating between different inffation models lies in fact that all cur- 
rent models of inflation predict the same near-scale invariant spectrum of scalar fluctuations. 
To further narrow down the number of observationally consistent inflation models, observ- 
ables independent of the scalar perturbation need to be measured and compared to model 
predictions. Two additional inflationary observables are the spectrum of tensor perturba- 
tions Pft, [21 E] and the non-gaussianity f^i [U El E] of the CMBR temperature perturbation 
spectrum. In recent years greater progress has been made in measuring these quantities 
directly. Upper bounds on the amplitude of the tensor perturbation spectrum, which is 
the spectrum of relic gravitational waves, have been determined directly through analysis of 
the CMBR polarization [71 18]. The non-gaussianity, which represents the deviation of the 
curvature perturbation from gaussian statistics, is also being better understood. Analysis of 
WMAP3 data [H] has found evidence of non-gaussian statistics in the CMBR temperature 
spectrum. With a better knowledge of these extra observables it becomes possible to better 
determine which model of inflation is most likely to have taken place. For example, a large 
Jml would tend to rule out a single field inflation model with minimal kinetic terms, while 
favoring those models that predict a large non-gaussianity. 

Ultimately, one would like to use the features of the CMB temperature anisotropy to 
reconstruct the inflaton action directly. It is customary to write the general scalar field 
action as [10] 



S = j d^x./^p{4>, X), X = ^g^^.d^d'^ 



Note that the Lagrangian (density) p is exactly equal to the pressure, which is the mo- 
tivation behind denoting the Lagrangian by p. In our analysis we will limit ourselves to 
actions that contain no third or higher derivatives of the inflaton. Throughout this pa- 
per we will assume that the curvature of the 3 non-compact space dimensions will be zero. 
Following the cosmological principle, our only choice for the metric is the FRW metric: 
g^" = diag(l, — a^, — a^, — a^), where a is a time dependent scale factor. Using ([l]), the 
Friedmann equations for the scale factor are 

?,MlH^ = p, (2) 

-2MlH = p + p, (3) 

where H = °^" is the Hubble parameter and p is the energy density, which in terms of the 
Lagrangian is 

p = 2Xp^x - p. (4) 

In single field inflation models with a minimal kinetic term the action is 

S= f d'x./^g[X-Vm. (5) 



If we assume (tsl), the only function that needs to be determined from the data is the potential 
V{(t)). Reconstruction of the inflationary potential for models of the form (|5| has been studied 



extensively pT| [T21 [131 [IH [ISl [ISl [IZ!- However, by assuming that the action has a minimal 
kinetic term we neglect a rich class of models such as DBI inflation [181 HH] , k-inflation [20] 
and ghost inflation [21]. In contrast, only a hand full of articles have been written that deal 
with the reconstruction of inflationary action with general kinetic terms [221 ESI l2l] . 

In non-minimal kinetic models the speed at which scalar fluctuations propagate can be dif- 
ferent than the speed of light. This can effect the temperature anisotropy in two ways. First, 
if Cs < c = 1, scalar fluctuations have a sound horizon that is smaller than the cosmological 
horizon, causing curvature perturbations to freeze in earlier than normal. Depending on how 
the Hubble parameter and sound speed change during the course of inflation, the temper- 
ature anisotropy can develop noticeable signatures of non-minimal kinetic terms. Second, 
models with non-minimal kinetic terms will in general produce a non-gaussian spectrum. 
Traditionally, the non-gaussianity is measured by the non-linearity parameter /jvl defined 
by the following ansatz for the curvature perturbation: 

c = Cl - hNiCi (6) 
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Here, ( is the general curvature perturbation and Cl is a curvature perturbation with gaussian 
statistics. Within the standard canonical action ([s]), non-gaussianities can be produced by 
cubic or higher order terms in the inflaton potential or by secondary interactions with gravity 
[5] . However, non-gaussianities produced by these mechanisms are on the order of the slow 
roll parameters, and thus small. In contrast, models with non-minimal kinetic terms can 
have large non-gaussianities, providing a clear distinction from canonical inflation. 

The goal of this paper will be to reconstruct an inflationary action from observables starting 
with as few initial assumptions as possible. In this paper we take the experimental inputs 
to be the scalar curvature perturbation Pg, the tensor curvature perturbation Ph, and the 
non-gaussianity (non-linearity) parameter /^vl. Unfortunately, completely reconstructing 
the off-shell action is not possible since the observables only carry on-shell information. 
To understand why the off-shell action is inaccessible to us, consider the interpretation of 
the action p{(j),X) as a surface in the three dimensional space (0, X, p) [23]. Because the 
observables are insensistive to the off-shell behavior of the action, we can only determine 
the one-dimensional trajectory p = p(0, X(0)) of the action on-shell, embedded in the two- 
dimensional surface deflned hj p = p{(j),X). A one-dimensional trajectory has an inflnite 
number of surfaces that contain it, each related to one another by a canonical transformation 
[23]. Therefore we have to be more speciflc about the form of the action that we are trying 
to flnd. In this paper we will reconstruct inflationary actions that have the form 

p{<P,X) = Pig,iX),...,gUX)JM,-Jnm. (7) 

Here it is assumed that P{xi, ...,Xm,yi, ■■■,yn) is a known function of the {xi} and {ya}, and 
the functions {gi} and {fa} are not all known. Once the on-shell trajectory = 0(/c) is 
determined, the action ^ deflnes a surface in the (0, X, p)-space up to a fleld redeflnition. 

The idea will be to use data on the CMBR perturbation spectrum to flnd the functions 
{gi{X)} and {fa{4>)}- In a naive comparison with algebraic linear equations, we expect 
that if there are n unknown functions, flnding the action requires n experimental inputs. 



Since we are assuming that there only three observables: Pg-, Ph and /atl, we can derive 
three reconstruction equations, which can determine an action with three or fewer unknown 
functions. In the case where the number of unknown functions is less than the number 
of experimental inputs, the reconstruction equations not used to find the action become 
constraint equations. 

Since we are interested in solving for functions {gi} and {/a} and not just numbers, we will 
need to know at least a portion of Pg, Ph and /at/Jj as functions of the scale k. While the 
scale dependence of Pg is known to be at least approximately power law dependent on fc, the 
scale dependence of the other two observables Ph and f^L is at this point unclear. Although 
future experiments will be able to clarify some aspects of the tensor and non-gaussianity 
signals, their exact functional forms will probably not be available for quite some time if 
at all. Regardless, the method we develop here does have utility outside of reconstruction. 
This method is well suited to testing how the form of an action depends on the observables. 
For instance if the scalar perturbation is of the near scale invariant variety: 

Ps oc P^'-i (8) 

we can use Pg to help derive an action and study its dependence on the index ng. That way 
if we wish to connect the action derived from (Is]) to an action derived from theory, we can see 
if the theoretical action leads to reasonable results for the observables. Furthermore, as we 
mentioned earlier when there are only one or two unknown functions in ([T]), the remaining 
reconstruction equations determine new consistency relations. In this paper most of the 
examples we deal with have only two unknown functions, which we solve for using the scalar 
and tensor spectrum data. The reconstruction equation derived from the non-gaussianity will 
then be a constraint; relating the non-gaussianity to the sound speed, the Hubble parameter 
and/or their derivatives. Outside of deriving the action we also find a method for quickly 
obtaining the sound speed as a function of time from the scalar and tensor perturbation 
spectra. Finding a sound speed different from the speed of light even over a small range of 
scales would be a powerful indication of non-canonical inflation. 

This paper is organized as follows. In section[2]we present the method for reconstructing the 
action from the scalar, tensor perturbation and the non-gaussianity parameter. We explain 
how cosmo logical data can find the Hubble parameter if, and the sound speed c^, and how 
these in turn can be used to find three unknown functions of the action ([t]). Once the 



method has been explained in section 2A_ we carry out a derivation of the action for different 
functions P(zi,Z2,^3), assuming that the scalar and tensor power spectra both scale like k 
to some power. In section |3] we apply our method to find the warp factor and potential in 
a generalized DBI inflation model. We find the warp factor and potential as functions of 
the spectral indices and the initial value of the Hubble parameter. The results for these are 
compared to the theoretically motivated warp factor and potential used in D3-brane DBI 
inflation. Finally, in section |4] we review our main results. 

^In this paper /jvl represents the equilateral bispectrum, and is therefore a function of a single scale. 



2 The Reconstruction Equations 

We start our derivation of the reconstruction procedure by explaining how the observables 
are used to find the Hubble parameter H and sound speed cj^ Once we have these, the 
action can be obtained using a set of reconstruction equations that will be shown later. Let 
us begin by recalling the definition of the slow roll parameteijj e in terms of the Hubble 
parameter. The definition of e implies that 

^ = -^H'- (9) 

dt ^ ' 

Since the perturbation spectra and non-gaussianity are functions of k and not time, we wish 

to rewrite this equation for ^ into an equation for H= ^^ff^- However, because we are 
assuming a general sound speed, we must be careful to differentiate between the horizons of 
scalar and tensor fiuctuations. If the sound speed differs from unity (in particular Cs < 1), 
then the horizon size of scalar fiuctuations: {aH/cs)~^ is smaller than that of the tensor 
fiuctuations: {aH)~^. This implies that at any given time, the scales ks and kt at which 
the scalar and tensor fiuctuations leave their respective horizons, will in general be different. 
For our purposes, we choose to study the dependence of H on the scalar wave number kg. 
Therefore, the condition for horizon exit is now kcg = aH instead of the more familiar 
relation: k = aH. Having made clear our choice of wave number, we now set out to express 
°^ in terms of familiar quantities: 



dt 



dt ' H dt 

where we have defined k = — — . Solving for °f above: 

Cs ^ dt 

dlogk _ H{l-e) 
dt 1 — K, 



dlogk ^^^_^^.d_logk^ ^^^^ 
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With equation (11) in hand, the equation for H can now be found from (9): 



He(l - k) 
H= ^3^. (12) 

We will use this equation to find Cg once we have found H and e in terms of the observ- 
ables. Since e depends on the time derivative of the Hubble parameter, H{k) and e{k) are 
independent parameters. Since we have two independent parameters, we will likely need two 
independent observables. The two observables we will use here are the scalar and tensor 
perturbation spectra. Recall that to first order in the slow roll parameters the perturbation 
spectra are given by 



Ps{k,) ^ 



87r2Mp2ec. 



(13) 

ksCs=aH 



^The method described here was inspired by the technique used in [2S] 

■^The term "slow roll parameter" is taken from chaotic inflation where inflation occurs only when the 
inflaton "velocity" (J) is small. However, DBI inflation can still occur for large </>. 



Phih 



2 H^ 



n- 



M. 



(14) 



pi 



kt=aH 



The extra source of information can also be garnered from the non-gaussianity parame- 
ter /iVL. However, going in this route would result in a more complicated solution. The 



parameter e can be found as a function of wave number using (13). Solving for e we obtain 



1 



H'' 



\ksCs=aH 



8nm^,Ps{k 



(15) 



ksCs=aH 



As a matter of convenience define Ps = AVs, where A is the value of the scalar perturbation 
at some fiducial scale k = k^. If ko ~ 0.002Mpc~^ then present data suggests that A ~ 10~^. 
Here, Vs is the normalized scalar perturbation defined such that Vsiko) = 1. Furthermore, 



let H = aH where a^ = Sn'^MpiA. Substituting (15) in for e in equation (12) we have 



n 



n' 



ksCs=aH 



c,,K - K' 




(16) 



ksCs=aH 



We have eliminated e from (12), but two independent variables remain. To get an equation 



for Cs we need to find 7i in terms of the observables. Since the expression for Ph ( 14 ) only 



depends on H, it can be used to find the Hubble parameter directly. In doing so, one find 
that 



n' 



kt=aH 



Vhjkt 
16 



(17) 



where Vh = A Ph. This gives the Hubble parameter as a function of the tensor mode wave 
number kt. In order to find IH?' as a function of the scalar mode wave number kg, note that 
the relation between the wave number of tensor and scalar modes that exit the horizon at the 



same time is kt 



k^c^. Therefore, we can obtain 7i 



ksCs=aH 



n' 



VhiksCs 



\ksCs=aH 
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by performing the substitution 



;i8) 



Plugging this in for Ti into equation (16) we get 



dVh{kt) ^ 

dlogkt 16cs{ks)Vs{ks) - Vh{kt 



2Vh{kt 



(19) 



Rearranging (19) we find the equation for Cg'. 



Cg[ks 



1 VhjksCs 
16 Vsik,) 



2Vh{ksCs 
Vh (ksCs 



(20) 



Here, a solid circle over Vh will denote differentiation with respect to logkf, not log kg. Once 
we specify what Vs{ks) and Vh{kt) are, we can use the above to solve for Cg- Even if we 
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can only show that c^ 7^ 1, this will be a signal for non- minimal kinetic terms. Once we use 



equations (18) and (20) to get H{k) and Cs{k) we can use (11) to find the relation between 



the wave number and time and ultimately find Ti.{t) and Cs{t). 

Having successfully found Ti and Cg, the next step will be to use this information to find 
the action. The most general single scalar field action is a multivariable function of (p and 
X. However, since we only have Ti and Cg as functions of a single independent variable (in 
this case k) the action can only be determined as a function of k: p{k). To turn p{k) into 
p{(j),X), we need to find and X as functions of k, invert them to get k{(j)) and k{X), and 
substitute into p{k). However, there is an ambiguity in how we substitute k for and X. 
Whenever k appears in the expression for p{k), we do not know whether to substitute it with 
A;(0), k{X) or some combination of the two. The ambiguity can be partially resolved if at 
the onset we specify a separation of the action into functions that depend either entirely on 
(f) or entirely on X. In light of this fact we make an ansatz: 

p{(f),X) = a'M^ig{^,x) = a^M}iP{g,{x),...,gUx)Ji{v),-Jn{v)) (21) 

where ip = MZ^cp and x = {aMpi)~'^X. We have also defined a dimensionless action g in 
order to keep the exposition neat and clear. Here it is assumed that P{yi, ■■■,ym, zi, ..., Zn) is 
a known function of the {yi} and {za}, and the functions {gi} and {fa} are not all known. 



We say that (21) only partially resolves the ambiguity, because it is possible in certain 



circumstances that a field redefinition can leave the form of (21) unalterated. For example 
consider the action 

g{^,x) = f{v)g{x). (22) 

Under a general field redefinition (f = h{(p) the action 

Q{v,x) = f{hmg{{h\^)fx). (23) 

If the function g is such that g{x ■ y) = g{x) ■ g{y) then 

g{^, X) = f{hmg{{h\0)y)g{x) = f{^)g{x). (24) 

Thus, not all choices for the function P lead to unique separation between the functions of 

X and if. 

The first equation for the reconstructed action will be obtained from the definition of the 
sound speed, which is given by 

^Xp^xx + P,x 



Assuming that the action has the form (21), equation (25) can now be used to find a dif- 
ferential equation for the gis as a function of time. After some work, this equation is given 
by 



{gP) X /2xx 1 
2 X \ x"^ cl 



m='-^-A-i^ + i,-^]-i9P9) (26) 



where a dot denotes different at ion with respect to the dimensionless time r = at. Here, we 
have used the short hand notation: 



{gP) = J29^P^^ (^P) = 5^^,P,, 



(gPg) 



i=l 



i=l 



«J=1 



U' 



(27) 



where P, and P,- are P 



V 



g and P., 



dgtdgj 



Equation (26), however, is incomplete since x 



is not known exphcitly. To turn (26) into a more usable form, we need to obtain an equation 
for X and its derivatives in terms of the known quantities e and Ti. To find such a formula 
lets write out the expression for the energy density of a general single scalar field action ^ : 



p = 2xp^^ - p. 



dH . 



As a consequence of the Friedmann equations, p + p is proportional to ^^ 

dH p + p 



dt 



2M2 



(28) 



(29) 



Solving for p + p in (|28|), and substituting the result into (|29|) we find that 



Ml dH 



X 



eMl,H' 



P,x dt 



P. 



(30) 



where in the last step we have used the definition of the slow roll parameter. With (21) as 
our assumed form of the action, this algebraic equation for x becomes a first order differential 
equation: 



X 
X 



ehe_ 

{gpy 



(31) 



Differentiating this equation, one can find an expression for x. After substituting the results 
of these relations, equation (26) becomes 

(gP) 



2ec^n^ 



1 + c 



/ 1 " 



(32) 



where t] 



Ht 



and (gP) ^ denotes partial differentiation of the quantity (gP) with respect 



to fa- We refer to (32) as the sound speed reconstruction equation. The non-gaussianity 
parameter f^i can also be used to find an equation for the functions Qi and fa- Following 
from the ansatz of the curvature perturbation (|6]), /nl is determined by the behavior of the 
curvature three point function: 



(C(kOC(k2)C(k3)) = -i2nY5^'\k, + k, + ^s)PsiKf-^^^ 



(33) 



where ki = |kj| and K = ki + k2 + k^. As we can see from (33), the /nl "will depend on the 
size and shape of the triangle formed by the three scales of the three point function. In [6] 



the authors found an expression for f^L for general single field actions ([I]) when the three 
scales form an equilateral triangle: 



/ 



NL 



where 



35 

108 



A 



i_A-Aa-l-2A 
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ci 



Xp^x + 2X^p,xx 



(34) 



(35) 



To get fxL as a function of K, we have to evaluate (34) at the time when the scale K 
passes outside of the sound horizon: Kcs = aH. The scalar power spectrum depends on a 
single scale k, which has a one-to-one mapping with the time through the relation kcs = aH. 
However, since f^L really depends on three different scales, the mapping between time and 
scale is not as straight forward. When the delta function in (33) is taken into account. 



the non-gaussianity still depends on three degrees of freedom: the magnitude of two of the 
scales and the angle between them [26]. To simplify matters, two of these three degrees 
of freedom will be fixed, so as to make fxL a univariate function. Since the equilateral 
configuration: ki = k2 = k^, has been very widely studied O |26], we will take f^L to be 
the non-gaussianity of the equilateral bispectrum. The equilateral non-gaussianity will be 
a function of /catl; which is the length of the sides of the equilateral triangle. Since the 
non-gaussianity freezes in when the scale K leaves the horizon, the scales at which the non- 
gaussianity and the scalar perturbation freeze in are not the same but instead related by 
3kNL = kg. After some work, one can use (34) and (31) to show that the g^s and /'s satisfy 
the equation 



(gP) 



WeH^ct 



1 



c;c; T _ r2 — sIAr'if 



INL 



X \ K + 



{gPf 



Yl f- \(9P) ( {-9P),a + iaPg),^ ) - {gP),. ( CgP) + (gPg) 



(36) 



a=l 



where k 



Cs 

Hcs 



1-K ■ 



This is the non-gaussianity reconstruction equation. Note that 
(|36|) is only well defined if the last line is nonzero. If the last line does vanish and the right 
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275 s 



CsInl = 0, and the non-gaussianity (34) only 



hand side of (p2|) is nonzero then 
depends on the functions gi and fa through the sound speed Cg. We will discuss such a case 
in the next section. Finally, another relation between the /'s and g^s can be derived by 
combining the Friedmann equations ^ and (pi): 



,gm,fi,...,fn) = {2e-3)n' 



(37) 



The upshot is that we now have four equations: (31), (32), (36), and (37), which when 



combined can be used to find x{t) (and by extension (p{t)) and three of the functions fa 
and gi. With more observational inputs it may be possible to determine even more / and 
g functions, but for now we will be content with what we have. In what follows, we will 
consider different, specific scenarios for the action and show how the action in each can be 
determined from the observables. 
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2.1 Examples 



In the case where the action (21) has only one function of x the equations (36) and (32) take 
on a much simpler form: 



2ec2?/2 



'-Tfp) 

9 a=l ) 



{3t 



gPg 



IGeH^ct 



ecln 



55 



972 
275 ""s 



clf. 
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-+^E/^[^.^^ 



p^ 



gga 



PgaJ^^gg\ 



(39) 



a=l 



As we mentioned earlier not all forms of the action will yield an equation for the functions g^^ 
and fa- In particular if the action is such that PgPgga 



PgaPgg for each a, and if the sound 



speed is constant, then (39) is not well defined. To see why, lets assume that Cg is constant. 



As a result, according to the definition of the sound speed: 



Q.- 



^•^Q,xx ~r Q^x 



Q{^,x) = fi{ip)x^^ + f2i^), 



(40) 



where the /i and /2 are integration constants and in general will be functions of ip only. We 
already know that with this form of the action PgPgga = PgaPgg- Thus the term in (39) 



inside the large parentheses vanishes, however the right hand side of equation (38) does not 
vanish. Therefore, we expect the denominator 1 — c^ — ^c^fNi in equation (39) to vanish. 



Indeed if we use the formulas ( 34 ) and ( 35 ) for f^i we find that 



I 



275 /I 



NL 



972 V c" 



(41) 



This relation between Cg and /nl holds regardless of what the functions /i and /2 in (40) 
are. It might be argued that if Cg is constant then the action (40) can be assumed and 
the remaining equations can be used to find /i and /2. This however is not that case since 
we have already used the sound speed equation to find g{x). This can be confirmed if one 



assumes the action (40). With (40) as our action, equation (32) is equivalent to the time 



derivative of equation (31). Thus, there are really only two equations: either (31) or (32) 



and equation (37). Therefore, only one of the two /i and /2 can be solved for. 



There is still yet another potential complication that may arise, specifically when the action 
takes the form g{(p,x) = f{(p)g{x)- Using (37) to find / in terms of g, the equations (p8J) 



and (39) become 



^n-/3„2 



GerjH c 



gf 



^^ 3(l + c2)-2e' 



55 



_ 972 2 f 

275 ^s J ' 



NL 



(42) 
(43) 
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As with the previous case, (43) is not defined when k = 0. Furthermore, the first equation 
(42) is also undefined when fj = 0. Since this is equivalent to e = constant, lets assume that 



e is constant to see which type of action this corresponds to. From (30) 



P 



xg , 



Q 



Q{ip,x) =/i((^)a;2-3. 



(44) 



Notice, that this action is a special case of the c^ = constant action (40) with /2 = 0. 
Thus, e = constant implies that Cg = constant. However, the converse of this is not true if 



/2 7^ 0. Since (38) is a well defined equation even with the action (44), we suspect that the 



denominator in (42) vanishes. After calculating the sound speed with the action (44) one 
finds that 

cl = -^^- (45) 



So indeed, the denominator in (42) does vanish. Assuming that the equations (42) and (43) 



are well defined, consistency requires that the right hand sides of these equations be equal, 
leading to the relation 



Sry 



in 



1 



3(l + c2)-2e 55 1 



972 
275 



clf. 



(46) 



NL 



This is a consistency relation between f^L-, Cs and the slow roll parameters. Although 
this consistency relation only holds for models with the action q = f{(f)g{x), analogous 
consistency relations can be found for any model in question. In what follows, we will carry 
out the full derivation of the g and / functions for two special cases. 



2.1.1 Case 1: g{ip,x) = g{x) -V (if) 



Suppose the action has the form 



q{'^,x) = g{x) -V{ip). 



(47) 



This type of action corresponds to the standard scalar field action when g is the identity 
map: g{x) = x. We refer to g{x) as the kinetic function. Notice that we have replaced what 
should be /i in our previous nomenclature with V{ip) in order to draw a clear analogy with 



the potential in the canonical scalar field action. With this type of action the equations (38) 



and (39) become: 



9 



2eclre (fj - 2e) 



9 



1 + C2 

IGekH^ct 1 



55 



1 _ ^2 - ^Jlr'^f^rr ' 



(48) 
(49) 



Here we have two expressions for the derivative of 5'(t). Consistency demands that the right 
hand sides of these equations be equal, thus we are lead to an analogue of the relation ( [46| : 

"-^'- *** ^ (50) 



l + c2 551-d-i|d/, 
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Interestingly enough this is the same consistency relation found in [23] for general single field 
infiation models. However, while the relation in [23] was only approximate, in our case it is 
exact. This shows that the consistency relation of Bean et al. is exact in the case when the 



action is of the form (47). Continuing with the derivation, the equations for V{t) and x{t 

;ctiv 



are given by (37) and (31), respectively. They read 

g{r) - (2e - 3)7^^ 



X 

Hx 



eW 



l + c' 



(51) 
(52) 



With our equations in hand we are almost ready to solve them and find the action. However, 
we still lack knowledge about the Ti and c^. In order to go further we need to look back to 



section^ and in particular equations (18) and (20). In order for these equations to be of any 



use we need two observables as inputs. For these we will assume that the two inputs are the 
scalar and tensor contributions to the CMBR. Presently, it is believed that these spectra are 
near scale invariant, and over a limited range of scales possess the formg^ 



Vsik) 



,{ns--l)logk/ko 



'p.(k) = Q^(nT-l)\ogkcs(k)/koCsO 



(53) 
(54) 



Here Cso = Cs(A;o) where fco is the fiducial scale at which Vs = 1. Note that we have assumed 
that the spectral indices have no running: i.e. n^, ut = constant. Although running spectral 
indices is an interesting generalization, in order to better demonstrate the usefulness of this 



procedure we will stick with simpler case of no running. Equation (18) then tells us that the 



Hubble parameter is simply proportional to the square root of the tensor perturbation: 



n{k) 



nj_ l^gf^C^^^l^yj^g^^g 



(55) 



4 4 

Defining Hq as 7i(fco) = T^o, the constant B is therefore B = IGTYq. As it stands, (55) is not 
complete since we still do not have an expression for Cs{k). To find c^ we turn to equation 



(20), the solution of which gives us 



CsO 



Cso nT-l\"T 



log fc/fco 



Since we defined Cgo as Cs{ko 



Hlnx — 3^ 
CsO, consistency of our definition demands that 



CsO 



^2^T-3 
rtn — 



(56) 



(57) 



nx — I 

Note that if 1 < n-r < 3, Cso is negative: a nonsense result. Therefore, we must restrict n^ 
to be either n^ < 1 or n^ > 3. With an expression for Cs in hand, 7i{k) explicitly in terms 
of k is 



n{k) = Hoe 



^"v^L":--| i°gfcAo 



(58) 



'*Note that the scale invariant tensor spectrum is most commonly defined at ut = 0, whereas we have 
defined it at ut = 1. The reason for using a nonstandard definition of ut is to establish a parity between the 
forms of the scalar and tensor perturbation spectra. The relation between the two different n^'s is simply 



riT 



standard 



"-T 



1. 
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Note that e and k are constant in this case and are given by 

re _nT-l 
CsVs nr - 3 ' 



(59) 



c^ riT — n^ , , 

K= — - = — -. 60 

Cs nT-2 ^ ' 

Since e is a constant then r/ = 0, which will simplify matters later when we try solve the 



reconstruction equations. Solving for log fc in ( 11 ), we find that the time dependence of log k 
is 

log k/k^ = - log [1 + m,{r - To)] , (61) 

K 



K 1-e _ 2(nr-ns) 

e 1 — K 

functions of time are 



where uj = f jz^ = — / _f-)("_2) - Therefore, the sound speed and Hubble parameter as 



1 + e/to(r - To) 

These are the expressions for the sound speed and Hubble parameter that will be used 
throughout this paper. They are completely independent of the form of the action that we 
are solving for, and are determined only by the inflationary observables Vs and Vh- 

Before we go about solving the reconstruction equations we should point out that not all 
values of the spectral indices lead to realistic inflationary scenarios. As has been mention 
before, the sound horizon of the scalar fluctuations is not the same as the cosmological 
horizon. As a consequence it is now possible for the size of the sound horizon to increase 
as time progresses. Thus, the usual expectation that larger scales freeze in at the beginning 
of inflation and smaller scales freeze in at the end, is not always guaranteed to hold. Recall 



that the time dependence of the scale is given by equation (61). It follows that the sound 
horizon size depends on time like 

Sound Horizon Size oc {-yaH)'^ = {k/ko)~^ = (1 + eHo^r — tq))~^ . (63) 

If uj/k, > 0, the sound horizon decreases with time as is normally expected. However, if 
uj/k, < 0, the size of the sound horizon increases during inflation, allowing modes the were 
previously frozen-in behind the horizon to reenter while inflation is still going on. This 
would be a disaster since if the scales were to reenter during inflation they would continue 
to fluctuate, destroying the near-scale invariance of the CMB anisotropy. Therefore, the 
spectral indices must be fixed such that uj and k, are either both positive or both negative. 

Since we are considering only those models that allow for inflation, we need to be sure 
that the spectral indices are such that an inflationary phase is allowed. If we refer to the 
expression for the equation of state w we see that not all values of ut are allowed if we want 
to have inflation: 

P P nT-7 

p 2Xp^x - P 3(raT - 3) 
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Notice that so long as nj^ < 3 the equation of state is always w < —^, and so inflation will 
occur. Since Cg oc e~\ in order for q to be interpreted as a sound speed, e must be positive. 



If we look back to equation (59) we find that not all values of ut will result in a positive 



value for e. Requiring that e > 0, we find that ut must be either n^ < 1 or n^ > 3. Since 
we have already found that ra^- > 3 would not lead to an inflationary solution, we conclude 
that riT < 1- Recall that in the previous paragraph we found that the sound horizon could 
expand during inflation only if the spectral indices were chosen so that uj/k, > 0. If one refers 
back to the definitions of u and k in terms of the spectral indices, we can see that if nr < 1 
the scalar spectral index is required to be Ug < 2. 



The sound speed (62) can tell us something about the expected range of validity of the 



scalar (53) and tensor spectra (54). If a; < then at some time t > tq the sound speed will 



be greater than one, signaling that fluctuations propagate at superluminal speeds. Likewise, 



superluminal speeds also occur at times t < tq when u > 0. Therefore, (53) and (54) can 



only be considered approximations; reliable within a certain range of scales. Keeping in mind 
that the sound horizon needs to shrink during inflation, the wave number k must respect the 
following bounds if the sound speed is to be less then the speed of light: 



k/ko > (c,o) I"! ^ 

i_ 

k/ko < (cso) i«i . 



for uj > 0, 
for uj < 0. 



(65) 



The only way (53) and (54) could be acceptable at all scales is if cu = 0, in which case Cg 



is a constant. If it turns out that the sound speed is not constant, (53) and (54) are most 



likely too naive. The most recent data from WMAP [S] suggests that the scalar spectral 
index may have a small but nonzero running, so we should not be surprised that our simple 
expressions for the perturbation spectra are not exactly correct. Regardless, scalar and tensor 
spectra with constant spectral indices are still a good approximation to the CMBR data. 
Our discussion will still be of relevance, as long as we keep in mind that the reconstructed 
actions are only approximations, valid over a limited range of scales. 

We will now simplify our discussion by fixing the sound speed to a constant, which is 
achieved by setting u = 0. Although we will be considering only constant sound speeds, we 
will keep the value of Cs arbitrary. This will allow us to find a more general solution to the 
reconstruction equations, while allowing us to study the limit c^ ^ 1 and see whether the 
canonical action is recovered. One might object to this choice of u on the grounds that if 



Cs is constant the non-gaussianity reconstruction equation (49) will be ill-defined for reasons 



discussed in section |2.1[ However, we counter that this is acceptable since we are assuming 
that only two functions g and V are unknown, thereby making the third reconstruction 



equation (49) unnecessary. It should be pointed out that while the two unknown functions 



can still be found when u = 0, the consistency relation (50) is no longer well defined. Once 



we substitute (62) for Cg and H into (48) and solve for g{T) the result is 



9{r) 






:i + enoiT-To)y 



90- 



(66) 
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Using the Friedmann equation (51) we can now find V{t): 



V{r) 



TLl 



3(1 + 0^0) -26 ^^^2 



Similarly, the equation for x is given by 



2ec,o +^0- 



(67) 



X 



Aec 



sO 



Tin 



X l + c%[l + eno{T-ro)]' 



(68) 



and the exact solution for x is 



=?o 



x{t) = Xo [1 + eHir - tq)] ^+=^0 



(69) 



Integrating (69) to find (p{t) 

Vo l + c^o 
eUo 1 - c2 



^ 



sO 



:i + e7^o(r-ro))^+=^o -1 



+ </2o, 



(70) 



where cpo = V^^o- Now that we have x and cp as functions of time, we can invert these and 



substitute the results into (66) and (67) to find g{x) and V{(p). After carrying this out, we 



can combine the g{x) and V{(p) to arrive at the full action: 



g{'^,x) 



n nj2„2 



l + '^sO 



nl[3{l + cl) - 2e] 



1 



-sO 



l + eH. 



CsO^ 



Vo 



1 + C2 



sO VO 



,2(1+4)1 



(71) 



Here is the complete action in the case when the sound speed is constant. Note that the 
final result does not depend on the integration constant qq. This is a result of the fact that 



the right hand side of equation (37) is independent of the initial values of the kinetic and 



potential functions. The only undetermined constants are the initial values of the scalar 
field and it's derivative, and due to the attractor nature of infiation, their exact values are 
unimportant. Despite what was said earlier in regards to the indefiniteness of the non- 



gaussianity reconstruction equation (49), this action does have a definite non-gaussianity 



given by the result in equation (41). It is worth noting that in the exceptional case where 
c, -^ 1: 



g{ip,x) = enl{x/xo) -nl{3-e)e *o 



^(v-vo) 



(72) 



we recover the canonical infiation action with an exponential potential. If we require that 
the wave function retain the standard normalization then xq = ehC^ and the action becomes 



q{^,x) 



X 



nlis 



e e 



~V2e((p-(fo) 



(73) 



which is the action of power law infiation [27j. This is a reassuring result; it confirms that 
in the appropriate limit, we can recover the standard infiationary action. 
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2.1.2 Case 2: q{^, x) = f{v)g{x) - Vi^p) 



Let us now take the complexity of the action one step further and assume that there are now 
three unknown functions: g, f and V. We define g as 



Qi^,x) = f{^)g{x) -V{^). 



(74) 



In equation (38) the only nonzero Pga is the one corresponding to /. Thus (38) reduces to 



gf 



(75) 



Furthermore, with this action the terms with the /^'s in (39) all vanish. The final result is 
simply: 

IQeRH^cl 1 



gf 



55 



972^2 



275 ^s 



clf. 



NL 



Combining equations (75) and (76), g decouples and we get equation just for /: 

_^ ^ ^ ^~^ l + cl 



V 



2e 



55 1 - c2 



972 „2f ■ 



(76) 



(77) 



Once we have solved for / here we can substitute the solution into equation ( 76 ) and solve 



for g. With the solutions for these two, V is found using the Friedmann equation (37). The 



final step is to find (p{t) and x{t) by solving (31): 

. IQkHct 



X 



X 



55 1 



972^2 



"s 275 *"« J^L 



{7i 



Let's assume that k 7^ 0, so that the reconstruction equations (76) (77) and (78) are well 



defined. We will again assume that the scalar and tensor perturbation spectra are given 



by (53) and (54). Therefore, H and Cg are the same as those that we found earlier (62). 



However, now that we are using the non-gaussianity reconstruction equation we need to 
specify /nl- In this example we will take /nl = to simplify the analysis. With these as 
our inputs, the reconstruction equations become 



Hf 



-2e 



Secu 1 + ci 



.2' 



gf 



IQe^u c^,n' 



55 1 



.2' 



X 

Tix 



IQeuj 



55 1 



^2' 



(79) 



Here, we have used the fact that r^ = and k 
They are 



-tuj. Each of these has an analytic solution. 



f{r) = /o 



cl{r) 



-sO 



gij) = go + 



2nt{i-c-j)t. 



Csofo 



XiT 



uj 55 



c'Jt) 



1 - cljr] 



8/55 



^0) 



-sO 



Fic;\r)) - Fic;,') 



Xo 



.2 



c^r) 



(821 



16 




Figure 1: Plot depicting the function g{x). This plot was made with Cso = 0.5, e = 0.1, /o = 
1, 5^0 = and ipo = 1. The exact behavior of g{x) is contrasted against the approximation 



(85). The behavior of g{x) is very linear except for small deviations for x < Xq. Note that 



at a;o(l — c^qY^^^ ~ 0.48 the plot of the exact behavior of g{x) stops abruptly as a result of 
the fact that g becomes non-real in this region. 



Here we have defined F(x) as 



4 63 59 



where 2-^1 is a hypergeometric function. Note that we can find a complete expression for 



g{x). We simply have to solve for Cs{t) in (82) to get c^(x), which is 



clix) 



(1 



-55/8 



(84) 



and then substitute this for Cs{t) in (81) to get g{x). Interestingly enough, g{x) is indepen- 
dent of uj, so taking the cj — > limit here is trivial. Since the solution for g{x) is in terms 
of hypergeometric functions, to get a better idea of what g{x) looks like we expand around 
Cso = 1, and thus obtain 



9[x) = go + 



eHl X 



Xo 



enlil-cso) 



foCsO Xq 



2585/oc 



sO 



2145 - 2209— + 64 

Xo 



+ 0((1 - Co)' 



^5) 



Let's take a moment to comment on the analytic behavior of g{x). In fig. II] the exact 
functional behavior of g{x) is shown along with the approximate expression (85). As fig. 
flland the approximation (85) suggest, the behavior of g is nearly linear with respect to x, 
when CsQ is close to one. However, for reasons that will be clear shortly, the limit c^o —>■ 1 
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does not necessarily mean that the action will be linear in x. Another interesting feature 
of g{x) is that it becomes non-real for values of x less than a;o(l — c^o)^ • This implies a 



lower bound on the values of x, which is a behavior that is observed in the solution (82) 



This lower bound is a result of the fact that for x < xo(l — c^o)^''^^ ^^^ sound speed squared 



would be negative according to (84) 



As for the functions /(v^) and V{ip), one cannot find analytic expressions for these like we 
did for g{x). Once we integrate x{t) to find (^(t), we can see why: 



(p = (po + 



-sO 



N4/55 



niii + n: 



dir) 



^sO 



" 55 



W) 



-sO 



where we have again shortened things by defining F^ as 

4 59 



[X 



iFi 



,£ l_ 

'55 ^ 2uj'' 



55' 55 ^ 2cj' 



X] 



(86) 



^7) 



Since ip is such a complicated function there is no way to invert (86 ) to get time as an analytic 



function of ip. Therefore, we are forced to either evaluate /(v^) and V{ip) numerically, or 
make an approximation for T{ip). Since we will be interested in finding a correspondence 
with the example in the previous section, we will approximate v^(t) in the uj <^1 limit. The 
result of this approximation is 



r 



^0 + 



^sO 



bbeHc 1 



[(1 + eHor) log(l + eHor) - eHor] + 0{uj' 



^sO 



To get t((p) we will drop all uj dependent terms from (88), so that (p{t) is a linear function 



of r. This approximation turns out to be remarkably accurate even at late times, since the 



higher order terms in (88) scale only logarithmically with r. Now that we have at least an 



approximate expression for t (</)), /(v?) can be found by replacing Cs{t) with 

1 + ^{ip-ipo)] . (89) 

The exact behavior of /(v?) was evaluated numerically and the results are shown in fig. [2] 
Since we will be taking the u ^ limit later, we make a further approximation of f{f) by 
taylor expanding around a; = 0: 



fiv) 



/o 



:i + fi(^-^o))^ 



+ 



8^/ol + c^olog[l + ff(^ 



VO 



^o)] 



55 1 



;i + fi(v.-(^o))^ 



(90) 



Therefore, although the higher order terms 
— > 1 is taken the product fg will retain the 
nonlinear x terms. This is why the action may not be linear in x even when Cso = 1- To 



Notice that the second term diverges at Cgo 
in (85) vanish when Cgo 



1, when the limit c 



sO 



find the potential we use the Friedmann equation (37). Doing so requires us to find g as a 



function of (f, which we find by replacing Cs{t) in (81 ) with Cs{(p) (89). The potential, taylor 
expanded around u = 0, is 



V{V>) 



9ofo - (2e 



3)Hl 



(l + lf(^-^o))^ 



+ 



8uj go fa 



'l + c^ 



sO) 



2eH2c.olog[l + l^(^-V^o)] 



Lo'-sO 



VO 



55 



-sO 



(l + ff(v.-^o))^ 



(91) 



Vo 




Figure 2: Plot depicting the function /(v?) for different values of u. This plot was made with 
Cso = 0.5, e = 0.1, ifQ = 0, 00 = 1-0, /o = 1-0 and qq = 0. We have also included a plot of 



the function f{(p) in (95) for comparison. 



Note that the individual functions g, f and V depend on the arbitrary integration constants 
/o and go, even though the action that is composed of them does not. If we are interested 
in just finding the action, fixing /q and go would be a moot point. However, it does raise 
the matter of how one separates the action into kinetic and potential terms. For example, 
suppose we separate the kinetic function into a constant and a "x-dependent" piece: 



g{x) = c + G{x). 



(92) 



The constant c is arbitrary and can be adjusted to any given value by absorbing the difference 



into G{x). Substituting the right hand side of (92) for g{x), the action (74) becomes 

gi^, x) = f{v)G{x) - V{^) + c/(^). (93) 

With the action written in this way, it would make more sense to define G{x) as the kinetic 
function and define the potential as 



v{^) = V{ip) -cf{ip). 



(94) 



In the case where / is constant (such as the example in the previous section) then the above 
redefinition only amounts to a uniform shift in the potential. However, if / is non-constant 
then the behavior of the potential can change drastically. Although none of the CMBR 
data are sensitive to changes in c, it is possible to find a value for c by requiring that in 
the appropriate limit, the action becomes equivalent to the canonical action. We will define 
this as the canonical limit of the action. Before we determine c by this method we need to 



confirm that the action ( 74 ) is canonically equivalent to the canonical action ( 73 ) when the 



sound speed is constant and equal to one. 
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If we turn our attention back to our approximations for f{(f) and V{(f), we notice that 
taking Cso = 1 leads to divergent results. These divergence are understandable since the 



reconstruction equations (79) are divergent when c^ = 1. However, if we set a; = in (90) 



and (91), it's possible to take Cso = 1 and still obtain a well defined result. Doing so results 



in the following for the functions g, f and V: 

^(-)-^o + ^^^^, m-^—^ — ^^ ^(^) - Hx^^^^' (95) 

where now Hq = e since Cgo = 1. We refer the reader to fig. ^for a comparison of f{^) in 



(95) and /(v?) for general values of cu and Csq. The action in the c^o — ^ 1 limit when tu = is 

'^'^' ^ (i + ff(v.-v^o))^^o (i + f!^(v^-^o))^- ^''^ 



This action can be related to the canonical action ([73j) through the field redefinition defined 

by 



(y^-y^o) = eV 2(^-^0). (97) 



00 

Under this redefinition, the new action is 



This is the same as the canonical action (73) found in the first example. To ensure a smooth 



transition to the canonical action we must separate the kinetic function as we did in (92) so 



that G(0) = in the cj — > and Cso — > 1 limits. Upon inspection of g{x) in (95) we see that 
the redefined kinetic function G{x) is 

G{x) = six) - f^E^l^. (99) 

Jo 



In doing so the potential is redefined according to (94) as 



viv) = V{^) - ^°^°/^° /(y^). (100) 

Jo 

There is a subtlety in this analysis that should be addressed. In order to reclaim the canonical 
action we needed to take the limits tu — *> and Cso — ^ 1 simultaneously. In our case we took 
that limit by setting u = and then letting c^o approach one. However, this is by no 
means the only way to take the limit. For example, we could have approached the limit 
by setting u = 1 — c^q and then take the limit as c^o goes to one. Had we taken the limit 
from a different direction it is possible that the action that resulted could have been different 



from the canonical action (73). After some inspection, it can be shown that under a field 
redefinition ip = h{(p) such that f~^{ip) = {h'{(p)Y, the potential in the canonical limit is 
given by 

V{(p) = (3 - 6)H^e-'^(^-^«) + (^0/0 - en'^ ^. (101) 

Jo 
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Here /(<^) = f{h{(f)) is the function / when the canonical hniit is taken. It is simple to 
show that the canonical limit of / is not unique, which means that the potential is also not 



unique. However, if we redefine our potential according to (100) instead, the new potential 



V is unique, and the action that results is canonically equivalent to ( 73 ) 



2.1.3 Case 3: fi{^) and f2i'^) Unknown 

We now bring up a case that will be of particular interest to reconstructions of the DBI 
action. We start by assuming that the action g{(f, x) has the form 

^(<^,x)=P(/i(<^),/2(<^),x), (102) 

where /i and /2 are unknown functions of ip. Unlike the previous cases, the functional 
dependence of the action with respect to x is assumed to be known exactly. In this case it 
is possible to obtain a set of algebraic equations of the two unknowns /i and /2. The first of 
these equations can be most easily obtained by going back to the original definition of the 



sound speed (25): 



2 P,x , „ 1/1 



2 ^ ^^x p _ 1 p (103) 



This equation together with (31) and the Friedmann equation (37) are enough to find fi{^) 



and /2(v5) in terms of the observables. In the next section we will see explicitly how the 



equations (103), (31) and (37) come together to reconstruct the DBI action from the power 



spectrum data. 



3 DBI inflation 

In realistic string and M-theories, the number of space-time dimensions is 10 or 11 dimen- 
sions. The extra 6 or 7 dimensions are compactified to small sizes, leaving the effective theory 
at low energies a theory of physics in four dimensions. The various moduli that control the 
shape (complex structure moduli) and size (Kahler moduli) of the internal space, also deter- 
mine the nature of the four-dimensional low-energy effective theory. Therefore, fixing these 
moduli is an important step in establishing a connection between string theory and the stan- 
dard model. In recent years, much attention has been paid to flux compactifications as a 
potential means of stabilizing string modulQ In a flux compactification, various fluxes wrap 
around closed cycles in the internal manifold creating a potential for the complex structure 
moduli. The best known of the these takes place in type IIB string theory. Here the internal 
space is six-dimensional Calabi-Yau and the 3-form fluxes F^ and H^ create a superpotential 
that fixes the complex structure [29] . These 3-form fluxes source a warping of the geometry 
of the internal manifold. In the type IIB flux compactification, the ansatz of the line element 
is taken as 

dslo = h-'/\y)g,,dx^dx'' + h^'\y)gmndy'^dy^ . (104) 

'^For a review of flux compactifications see [28] , 
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Here h is the warp factor which is sourced by the fluxes and varies only along the dimensions 
of the internal manifold. In DBI inflation, which we will be considering in this section, the 
local geometry of the internal manifold is a Klebanov-Strassler throat geometry [30] , and is 
described by the metric 

g^ndy'^dy^ = dr'' + rHs\^. (105) 

Here ds\^ is the line element of a five- dimensional manifold X5, which forms the base of 
the KS throat. The coordinate r runs along depth of the throat. For our purposes we will 
only consider motion along r and integrate over the base manifold X^jj The warping of the 
internal space creates a natural realization of the Randall- Sundrum model |33], and has also 
provided model builders with a new approach to developing string theory based models of 
inflation. The most popular inflation model that makes use of this warping is DBI inflation, 
which is the primary focus of this section. 

In the simplest DBI inflation models [TSl EH ES] a -D3 brane travels along the r direction, 
either into or out of the KS throat. The D3 brane extends into the 3 non-compact space 
dimensions and is point like in the internal manifold. The standard DBI action for the 1^3 
brane is 



S, 



DBI 



Jd'xV^ [f-\<P)y/l-2f{<p)X- f-\<p) + V(0)] . (106) 



Here = y/Tsr (where T3 is the D3 brane tension) is a rescaling of the coordinate r and 
will play the role of the inflaton. The quantity f~^ = T^h~^ is the rescaled warp factor. The 



metric g that appears in ( 106 ) is the metric on the 3-1-1 dimensional non-compact subspace 
which describes the geometry of our familiar 4 dimensional space-time. We will continue to 
assume that the geometry of the 3-1-1 dimensional subspace is described by the FRW metric 
with zero curvature. The energy density and pressure in the non-compact subspace due to 
the brane are given by 

p = r' (7 - 1) + V, (107) 

p = iifr' (7 - 1) - V. (108) 

Here 7 is a new parameter, not found in the standard canonical inflation. In terms of the 
quantities in the DBI action, 7 is defined as 

7 = , '^ (109) 

Vl - 2/(0)X 

The 7 defined here is analogous to the Lorentz factor in special relativity, and will hence- 
forth be referred to as the Lorentz factor. The Lorentz factor places an upper limit on the 
speed of the brane as it travels through the KS throat. Since the kinetic energy of the brane 
is limited, this allows one to get a sufficient amount of infiation even with potentials that 
would be considered too steep to use in standard canonical infiation. 



^Fluctuations of the brane position along the transverse directions of the KS throat have been mentioned 
as a possible source of entropy perturbations |31|, 132] . These could serve as a further constraint on the form 
of the action. 
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In our study of the DBI model we will be assuming that the scalar and tensor spectra are 
approximately (53) and (54), respectively. With these as our inflationary observables, we 



found that e was constant (59). The fact that e is a constant indicates that inflation will not 
end on it's own, and instead some other mechanism such as D3-D3 annihilation ^36j must be 
used to provide a graceful exit. Since our study is concerned more with the physics during 
inflation, this topic will not be addressed further. We will now present a generalized DBI 
action, and show how it is reconstructed from the inflationary observables. 



3.1 A Generalized DBI Model 

Having sketched out the general method for reconstructing different types of inflationary 
actions in section [2| it is now time to apply these methods to a DBI-type action given by 



g{(fi,x) = P{x,J^{ip),V{ip)), 



(110) 



where 



P{Zi,Z2,Z3) 



-1 



(VI - 2^2^! - 1) 



Z3- 



:nv. 



Here J-'{^p) = a^M^J^ip) is the (dimensionless) warp factor in the throat, and V{ip) = 
{aMpi)~'^V{ip) is the (dimensionless) potential. In the KS throat geometry the warp factor 
is taken to be JF oc y^"^. The potential V is assumed by many to be quadratic in ip. For the 
purposes of this study we will not assume a priori any form for the functions JF and V, and 
instead allow the inflationary observables to determine them. Now that we have established 



the general form of the action, we can use the procedure outlined in section 2.1.3 to find JF 



and V. Turning to equations (103) and (31) we find the relations 



Pn 



T 



\-2Tx 



Solving for T and x: 



Hf^) 



- - 

2x Vc? 



2en^Cs 



X 



eH^Vl - 2Tx. 



(112) 



X = eH Cs- 



(113) 



Comparing the second equation above with (31 ) and recalling the definition of (109), we find 
that 



1 

1 



:ii4) 



This result is characteristic of DBI inflation and holds regardless of the warp factor and 



potential used. Having found T, equation (37) tells us what V is 



V(y?) = 37^2 






(115) 
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Having already found the expression for JF in (113), we can now write down the important 



reconstruction equations for V and T in terms of ?i, c^ and e: 



Tik) 



1 + c, 
2ecJH?' 



(116) 
(117) 



To turn JF and V into functions of ip we need to integrate our solution for x (113) to find ip. 



Taking equation (31) to find an expression for 0, we find that in the case of DBI infiation 



V 



118 



The sign of the right hand side of the equation is ambiguous, due to the square root taken 



to get this equation from (31). The sign is left arbitrary for now and will be specified later 



based on the requirement that if be positive. Once we solve for 9?(r) in (118) and invert to 
get t(v9), we can then find a solution for V{ip) and J^{ip)- Now that we have laid the ground 
work for generating the functions of the generalized DBI action, the next section will show 
how the perturbation spectra are used to obtain explicit expressions for J-'{(p) and V{ip). 



3.2 The Warp Factor and Potential in DBI Inflation 

In this section we will now use the program that was laid out at the end of the previous 



section to find an exact solution for the warp factor and potential in the action (110). We will 



again assume that the scalar and tensor spectra have a power-law dependence with respect 
to the scale k. Therefore, the sound speed, Hubble parameter and e are the same as those 



found in section 2.1.1 Thus, the potential as a function of time is 



V 



Hn 



2e 



(l + e7^o(r-ro))2 
and the warp factor as a function of time is 

^ = <^ + "-::. "" II 



n 



l + ^il + eHoir-To))- 



(119) 



2nt 



^ (I + eHoir - To)y 



■2uj 



(120) 



Since we want infiation to occur we need e to be such that e < 1, which according to equation 



(59) implies that n^" < 1. However, n^ = 1 can be dismissed as physically unreasonable since 



that would imply c^o oc - — i> oo. After substituting our solutions for 7i(r) and Cs,{j) into the 
equation of motion for '^[j') we get 



^ 



±72^2(1 + e7^o(r-ro))-t-i. 



(121) 



Once we integrate this expression we can obtain an answer for '^{t\ As a matter of conve- 
nience we will set the value of the integration constant that results to zero. Later, once we 
have found T and V, we will see that this choice allows for a correspondence between the 
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reconstructed functions and their theoretically derived counterparts. After integrating (121) 
we find that 



(^(r) 



Vdr = T'^-^^^{l + en,{r-r,)r-/\ 



:i22i 



eu 



Since we are interested in eventually connecting the reconstructed action with the standard 
DBI model we need to keep the infiaton, which is just a rescaled radial coordinate, positive. 
The sign that we choose in (122) will therefore depend on the sign of uj. We can write the 
general solution as 



^{'t) = Vo(l + e^o(r - To)) 



-1^/2 



where 



^0 



2^2 



Ho 



euo 



(123) 



(124) 



In the case where cu > the field ip decreases monotonically to zero as time passes, which 
implies that the brane is falling into the throat. This corresponds to the UV DBI scenario. 
If on the other cu < 0, then ip increases monotonically with time and the brane falls out of 
the throat, which corresponds to IR DBI infiation. Solving for time in terms of p 



1 + m^^T - To) 






(125) 



Substituting this for 1 + e7io{T — tq) in the expressions we found for the potential and warp 
factor we find that V as a function of p is 






2e 



\ 



]_ I '^0 / y 



V>() 



(126) 



and the warp factor as a function of p is 



^ 



^ 



2^0 V^o 



-2-- 



1 






(127) 



Furthermore, when 7 is expressed as a function of p, it takes on a very simple form: 

2 o 1 



7 = 7o 






9 9 ' 



(128) 



where 70 = — . It is interesting to note that (128) is the same as the approximate results 
found in the theoretically inspired DBI model |M[ . The potential and warp factor derived 
here are the same as those found in ^7\. There the authors reconstructed the potential 
and warp factor by assuming that the equation of state w = - was a constant and that 
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-a;/2 



ip (X T "^1^. In contrast, we have reconstructed the potential and warp factor under the 
assumption that the scalar and tensor perturbations are (53) and (54). The non-gaussianity 



in this DBI model is the same result that one comes across in the literature [6]: 



/ 



NL 



35 
108 



(129) 



This particularly simple result is a general feature of DBI, and is independent of the warp 



factor and potential. This result for the non-gaussianity (129) also follows from consistency 



of the reconstruction equations (32) and (36). Thus, (129) can be viewed as a consistency 



relation analogous to those found in (46) and (50). An interesting generalization to consider 
is 



Q{^,x) = P{g{x),T{^)y{^)) 



(130) 



(36). 



where P{zi, Z2, z^) is defined by (111). As a consistency check, one can easily show that if 



the non-gaussianity is equal to (129), and the potential and warp factor are given by (119) 



and (120), respectively, then g{x) = x is a solution to the reconstruction equations (32) and 



Having found the potential and warp factor as functions of the inflaton, we can now say that 
our task is at an end. Amazingly enough, despite the complicated form of the reconstruction 
equations an exact solution for V and JF was available even for semi-realistic scalar and 
tensor spectra. In the next section we will discuss the properties of the reconstructed action, 
and its correspondence with the theoretically derived action of DBI inflation. 



3.3 Discussion 



In section 2.1.1 we found that not all values of the spectral indices lead to inflationary 
and/or physically sensible actions. Specifically, we showed that unless ut < 1 the matter 
described by the action was unable to drive an inflationary phase. Furthermore, when this 
constraint on ht was considered in conjunction with the requirement that the sound horizon 
decrease as inflation occurs, the scalar spectral index had to be bounded like Ug < 2. These 
results hold for any reconstructed action that was derived assuming the observational inputs 



(53) and (54). However, even if these constraints are satisfied, it is not guaranteed that 



the reconstructed action is physically sensible when interpreted in the context of a given 
theoretical construction. For example, if we are to interpret the action reconstructed in this 



section as a DBI action, then 7 > 0. Doing so would be contrary to its definition (109) within 
the context of DBI inflation. In this case, since 7 = — the fact that we already have enforced 
the constraint Cs > automatically keeps 7 positive. We will see later, however, that the 



constraints found in 2.1.1 are not sufficient for our reconstructed action to be interpreted as 
a DBI action. 



First, let's consider the warp factor (127) and what constraints it places on the observables. 



Suppose ut < Tig. In this case the exponent of ip in front of the square brackets in (127) will 



always be negative. Therefore, for small values of ip the leading order behavior of T will go 
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Figure 3: Plot depicting the regions of the ng-riT parameter space such that J^ diverges at 
the origin. The region in the upper left hand corner is the region of IR DBI inflation while 
the corner at the bottom right hand corner corresponds to UV DBI. Only those values of Ug 
and riT with n^ < 2 and n-r < 1 were considered, since any point outside that region would 
lead to an unphysical and/or non-infiationary action. 



like if""" where a > 0. Thus, the warp factor increases as we fall into the throat, which is 
what we would expect for a warped compactification in string theory. On the other hand, 
if rtT > Ug the leading order behavior of JF will be (p to some positive or negative power, 
depending on the relative difference between Ug and tit- If the difference between ut and Ug 
is too small, then to leading order, J^ will scale like (p to some positive power. This indicates 
that the warp factor gets smaller as we reach the bottom of the throat, which is a scenario 
that is difficult to embed into a string theory compactification. However, If the difference 
between Ug and n-p is large enough, then it is possible to get a more sensible solution where 
JF ~ (y9~°. In general, the condition that JF increases as we approach the bottom of the 
throat implies that 



<0 



uj < —2 or cj > 0. 



'1311 



UJ 



The regions in the Ug-riT parameter space where the condition (131) is satisfied are shown 



in fig. |3| The region shown in light grey in fig. [3] is defined by u; > 0, or equivalently 
Ut < Ug. This region corresponds to the UV phase of DBI infiation. The region in dark 
grey is defined by cj < —2, or equivalently n-r > — ^;;^, and corresponds to the IR phase. 
As fig. |3] illustrates, if Ug is restricted to the presently favored value Ug ~ 0.96, then the 
value of Ut is tightly constraint in the IR region but relatively unrestricted in the UV phase. 
Furthermore, since ut < 1, only those models in the UV phase can have blue-tilt. Recent 
CMBR data favors a blue-tilted spectrum, but only if there is a running spectral index [8] . 
Although running spectral indices would be an interesting extension of this analysis, we will 
leave this topic to future studies. 



An unpleasant feature of the warp factor (127) is that it becomes negative when y9 > 
V^o'^o/v^- This is particularly distasteful since the metric (104) depends on /^/^, which 
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means that at sufficiently large (p the metric is imaginary. The values of r where the warp 
factor is positive are given by 



r-To > 



r-To < 






forcj>0 (UV), 
for cj < (IR). 



[132) 



One can show using equation (61) that this is equivalent to the bounds in (65). This is no 



coincidence; it is a result of the fact that T is proportional to 1 — c^ ( 113 ). Therefore, for the 



same reasons that were explained in section 2.1.1 



the perturbations (53) and (54) can only 



be used as approximations. The inequalities in (132) tell us that JF is a valid warp factor 
towards the end of inflation in the case of UV DBI, and at the beginning of inflation in IR 
DBI. The time tq at which the initial conditions are specified should be at the beginning 
of inflation in the UV scenario, and at the end in the case of IR DBI. If we choose Tq in 



this manner then the approximations for the perturbation spectra (53) and (54) will lead to 



realistic sound speeds for the entire duration of the inflationary episode. 



It is worth asking if the JF that we have derived in (127) can approximate the AdS warp 
factor derived from theory. It is clear from (127) that this can be achieved if and only if 
uj = 2. However, the only way we can get tu = 2 is if either i.) rix = 2 ii.) rij- —* oo 
or iii.) Ug = 1. As we have already seen, case i.) is unphysical, and case ii.) is difficult 
to imagine taking place. While case iii.) is unlikely to be true exactly, it is nevertheless 
the more realistic of the three, especially when you consider that observation suggests that 



0.96. If we do set Ug = 1, the warp factor and potential become 



H^) 



2 1 

7^ 



1 
2?' 



(133) 



V 



3e' 



-V 



l + ty.2- 



(134) 



In the case where the field range is smalQ these are approximately 

'■pi 1 



f\<P) 



2M2 



a2e4 d)^-' 



(135) 



V{ct>) 



e'a^'i - 2e) 



(136) 



Where we have reverted back to the standard, dimensionful /, V and (f) for clarity's sake. 
It is a bit of a surprise that it in process of trying to recover the AdS warp factor we have 
stumbled upon the commonly used potential in UV DBI inflation. If we take (135) and 



^In [38] it was shown that DBI inflation is only consistent when the magnitude of the inflaton fleld is 
sub-planckian: ip <^1. 
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demand that it is consistent with the theoretical result we can arrive at a condition on e in 
terms of the D3 charge. Recall that in the KS throat /(0) is given by 

m = ^, (137) 

where T3 = 77^3— t'tw and 



Consistency with (135) demands that 



2Mp2 _ nN 10^ 



a^e 



2e4 Vol(X5) m/^' 



(139) 



In order to get an inflationary phase A^ ~ 10^°, putting us well within the range of validity 
for the supergravity approximation. While it is interesting that the standard D3 brane DBI 
model can be recovered from a near scale invariant scalar power spectrum, it has been ac- 
knowledged that this inflation model is problematic. In [38j Baumann and McAllister found 
that while present bounds on non-gaussianity imply that A^ < 38, primordial perturbations 
imply that A^ > 10'^Vol(X5). These two limits are incompatible unless Vol(X5) < 10~^. It 
is not clear that such a space could be naturally embedded into a string theory compactifi- 
cation. More general warp factors and potentials have been considered in [22]. There it was 
found that models could not simultaneously satisfy bounds on the field range and observa- 
tional bounds on the non-gaussianity. Therefore, even though our warp factor and potential 
matches the theoretically based predictions, the problems inherent in the DBI model carry 
over into its generalizations. 



4 Conclusion 

In this paper we have presented a method for deriving the actions of single field inflation 
models using CMBR data. This method allows one to derive up to three unknown functions 
of the action using the scalar perturbation Pg, tensor perturbation Pf and the non-gaussianity 
/nl- After stating the reconstruction equations, we carried out the reconstruction procedure 
for two simple examples. For the purposes of the reconstruction, we assumed that the scalar 
and tensor spectra were power-law dependent on the scale k, with the spectral indices kept 
as free parameters. In the first example we assumed that the action had the form shown 



in equation (47), and used the reconstruction equations to obtain the action as a function 



of the spectral indices. In this example there were only two unknown functions, thus the 



reconstruction equations also lead to a consistency relation (50) between the /nl, Cg and the 
slow roll parameters. However, this consistency relation is only well defined when the sound 
speed is not a constant. 

In the second example, the action depended on three unknown functions and therefore 
required all three reconstruction equations. In order to simphfy the discussion we took 
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as our input for the non-gaussianity f^i = 0. Although we were unable to express the 
action in terms of elementary functions we were able to obtain the action numerically and 
approximately assuming Cso ~ 1 and a; ^ 0. We showed that the action in this example 
was canonically equivalent to the canonical action derived at end of the previous section. In 
discussing this example we also pointed out possible ambiguities in the program relating to 
how one defines a separation between the kinetic and potential terms. 

In section [3} we used the procedure to derive and study the warp factor and potential 
in a generalized DBI inflation model. Again, we assumed that both of the perturbation 
spectra scaled like k to some power. Exact expressions for the warp factor and potential 
were then derived, each having an explicit dependence on the spectral indices. The demand 
for a physically sensible DBI inflation model placed constraints on the spectral indices. In 
addition we found that the derived action approximates the original UV DBI inflation model 
in the case where -C Mpi. Unfortunately, the problems that have plagued UV DBI inflation 
are still present in our case. 

This procedure was shown to be useful in studying how the action of a general inflation 
model depends on the observables. For example, we found that if the scalar and tensor 
perturbation spectra went like k to an arbitrary power, the reconstruction would lead to a 
realistic inflationary model only if ny < 1 and Us < 2. Furthermore, to keep the speed of 
fluctuations from becoming superluminal, the range of k over which the approximations for 



the spectra (53) and (54) are taken, had to limited. When we reconstructed a generalized 
DBI action in section [3| further constraints were needed to keep the action compatible with 
an interpretation of DBI inflation. Specifically, we found that in order to keep the warp factor 
positive, the field range had to be limited. Furthermore, in the theoretically motivated DBI 
model, the warp factor increases as we reach the bottom of the warped throat. In order for 
this to be true in our reconstruction, the spectral indices needed to satisfy the additional 
constraints: uj < —2 or cj > 0. 

In this paper we have only considered the simplest of the DBI inflation models, which 
unfortunately suffers from several inconsistencies. However, there are many extensions of 
the D3 brane DBI model that can circumvent some of the problems of the original. Some 
of these extensions include using wrapped D5 branes [391 SO], multiple D3 branes [H], and 
multiple throats [19]. Each of these models has its potential advantages and drawbacks. 
Applying our reconstruction procedure may help to further elucidate their relative strengths 
and weaknesses. Furthermore, we have limited ourselves to perturbations with simple power 
law behavior. However, this naive assumption may be incorrect. It is easy to imagine that 
the spectral indices themselves are also scale dependent. Based on the results of this paper 
we can predict what kind of effect a running spectral index would have on the physics of the 
underlying models. For instance, in the generalized DBI model it is possible for the spectral 
indices to change during inflation in such a way as to pass from the IR to the UV phasqj 
Transition between phases would correspond to a completely different physical scenario, one 
where the brane falls out of one throat and back into another. Therefore, running spectral 
indices would describe multi-throat DBI inflation. A model which has so far been shown to 



^Fig. [s] implies that inflation can change between UV and IR phases only if it passes through the exactly 
scale-invariant point: irig^nj') = (1, 1). 
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be internally consistent 

This study has also raised some other questions that may be worth investigation. In 
particular what is the relation between actions that yield the same observables. It may be 
possible to define a group of transformations which leave the perturbation spectra and the 
non-gaussiantiy invariant. Such a set of transformations would allow us to classify actions 
based on the observables they yield. Another interesting possibility that came out of this 
study is the idea of using the reconstruction equations as a way to generate consistency 
relations between f^i, the sound speed Cg and the slow roll parameters. These questions 
will be left for future studies. 
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